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Abstract. 

Until now, HD 155448 has been known as a post-AGB star and listed as a quadruple system. In this paper, we study the system 
in depth and reveal that the B component itself is a binary and that the five stars HD 155448 A, Bl, B2, C, and D likely 
form a comoving stellar system. From a spectroscopic analysis we derive the spectral types and find that all components are B 
dwarfs (A: B1V, Bl: B6V, B2: B9V, C: B4Ve, D: B8V). Their stellar ages put them close to the ZAMS, and their distance is 
estimated to be ~2 kpc. Of particular interest is the C component, which shows strong hydrogen and forbidden emission lines 
at optical wavelengths. All emission lines are spatially extended in the eastern direction and appear to have a similar velocity 
shift, except for the [Oi] line. In the IR images, we see an arc-like shape to the northeast of HD 155448 C. From the optical 
up to 10 fim, most circumstellar emission is located at distances between ~1.0" and 3.0" from HD 155448 C, while in the 
Q band the arc-like structure appears to be in contact with HD 155448 C. The Spitzer and VLT/VISIR mid-IR spectra show 
that the circumstellar material closest to the star consists of silicates, while polycyclic aromatic hydrocarbons (PAH) dominate 
the emission at distances >1" with bands at 8.6, 11.3, and 12.7 /mi. We consider several scenarios to explain the unusual, 
asymmetric, arc-shaped geometry of the circumstellar matter. The most likely explanation is an outflow colliding with remnant 
matter from the star formation process. 
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1. Introduction 

HD 155448 (alias SAO 208540, or HIP 84228) was described 
as a B9 object in the "Catalogue of two-dimensional spectral 
types for the HD stars, Vol. 3" (Houk 1982). Three other mem- 
bers of this system are known according to the "Catalogue 
, of the components of double and multiple stars (CCDM)" 
(Dommanget & Nys 1994). The system's Hipparcos parallax 
of 1 .65 mas (Perryman et al. 1 19971 1 would correspond to a dis- 
tance of roughly d = 600 pc, but the uncertainty of the parallax 
is large (1.99 mas). 

Up to now, all studies have described this system as old and 
evolved. In their sample of 42 IRAS point source objects, 
Van der Veen et al. ( fT989l considered IRAS 17097-3210 
(alias HD 155448) - based on its IRAS colours - as a tran- 
sition object from the AGB to the planetary nebula (PN) 
stage. These authors also speculated about PAH emission at 
7.8 yi/m. However, in a survey of 2703 IRAS point sources at 



Dwingeloo, Effelsberg, and Parkes, examining the 1612 MHz 
OH-transition - which is one of the strongest lines for AGB 
stars, but decreases towards the PN stage - Te Lintel Hekkert 
et al. (119911 1 did not detect emission from HD 155448. 



Send offprint requests to: oschuetz@gmail.com 

* Based on observations collected at the European Southern Obser- 
vatory (programme IDs 075.C-0091(B) and 075.C-0091(C), together 
with archival and technical data) 



Malfait et al. (1998) points out that, in order to distinguish 
between Herbig Ae/Be stars and post-AGB stars, one should 
also consider the ([B - V],[U - B])-colour diagram, since post- 
AGB stars show an important Balmer discontinuity, while the 
regions populated by both classes overlap in infrared colour- 
colour diagrams. They also show that several other currently 
well-accepted Herbig Ae/Be stars had previously been classi- 
fied as transition objects between the AGB and the planetary- 
nebula stage. From the absence of a prominent near-IR excess 
(the excess only starts beyond 4.5 yt/m), these authors consider 
HD 155448 to be a Vega-type candidate. Spectral energy dis- 
tributions for the integrated fluxes have been given by Van der 
Veen et al. ( fT989l and Malfait et al. (fl99"8l 

Several important issues remain unclear: 

- What is the evolutionary status of this system? Does the ex- 
cess emission point towards a pre-main sequence stage in 
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Table 1. Summary of all imaging data presented in this work. 



Passband 


Instrument 


Date 


Pixel scale 




Airmass 






of nf crnt Q hpcnn 

Ul 111 LI. Ill S Ut.fi ill 


L J 


[sec] 




B 


EFOSC2 


Feb 26, 2006 


0.157 


8 


1.3 


V 


EFOSC2 


Feb 26, 2006 


0.157 


8 


1.2 


R 


EFOSC2 


Feb 26, 2006 


0.157 


8 


1.2 


J 


SOFI 


May 07, 2004 


0.144 


100 


1.0 


H 


ADONIS 


Jun 08, 2000 


0.100 


600 


1.0 


H 


SOFI 


May 07, 2004 


0.144 


100 


1.0 


NB 1.64 fim 


NACO 


Sep 09, 2005 


0.027 


420 


1.1-1.5 


SK 


ADONIS 


Jun 08, 2000 


0.100 


600 


1.0 


Ks 


SOFI 


May 07, 2004 


0.144 


100 


1.0 


NB 2.12 jum 


NACO 


Sep 09, 2005 


0.027 


435 


1.1-1.5 


PAH2 (11.25 pm) 


VISIR 


Apr 27 - Aug 05, 2005 


0.075 


5425 


1.0-1.5 


PAH2_ref (11.88 /mi) 


VISIR 


Apr 27 - Aug 05, 2005 


0.075 


5425 


1.0-1.6 


Q2 (18.72 pm) 


VISIR 


Aug 26 -Sep 01, 2005 


0.075 


7948 


1.0-1.2 



Table 2. Summary of all spectroscopic data presented in this work. See Fig. Q] and Appendix IA.2.21 for the definition of the 
EFOSC2 slit orientations. 



Passband 




Instrument 


Date 


Slit 


Spec. res. 


*exp 


Airmass 








of night's begin 


orientation 


R 


[sec] 




Grism 5 


(5200-9350 A) 


EFOSC2 


Aug 19, 2007 


East- West 


-300 


1700 


1.3 


Grism 1 1 


(3380-7520 A) 


EFOSC2 


Mar 01, 2006 


A-C, B-C, C-D 


-400 


240 


1.2 


Grism 1 1 


(3380-7520 A) 


EFOSC2 


Apr 20, 2006 


East- West 


-400 


2700 


1.0 


Grism 18 


(4700-6770 A) 


EFOSC2 


Aug 19, 2007 


East- West 


-600 


2000 


1.1 


Grism 19 


(4441-5 114 A) 


EFOSC2 


Aug 12, 2008 


A-C, B-C, C-D 


-3000 


420-1800 


1.0-2.0 


Grism 20 (6047-7147 A) 


EFOSC2 


Aug 12, 2008 


A-C, East- West 


-2500 


420-1800 


1.3-1.8 


8-13 ytim 




VISIR 


Apr 15 - Sep 03, 2005 


East- West 


-350 


1800 


1.0-1.4 


SL 


(5.5-14 jum) 


Spitzer/IRS 


March 22, 2005 


see Fig.[3] 


-160 


12 


n/a 


SH 


(10-20 pm) 


Spitzer/IRS 


March 22, 2005 


see Fig.[3] 


-160 


60 


n/a 


LH 


(19-38 /an) 


Spitzer/IRS 


March 22, 2005 


see Fig.f3] 


-160 


24 


n/a 



which the star is surrounded by a disk? What is the compo- 
sition of the material causing the IR excess emission? 

- Up to now it has not been proven whether the components 
are physically related. Are there more members to be dis- 
covered? 

- What is the actual distance of this system? 

In this paper, we present the first high-resolution, multi- 
wavelength photometry and spectroscopy of the HD 155448 
system intended to address the above questions. In Sect. [2] we 
derive new spectral types for all components and constrain the 



system's evolutionary status and distance. In Sect. [3] we show 
that the B component itself is a binary, uncover arc-shaped IR 
emission near the C component, and derive the composition of 
the material emitting in the IR. After a discussion in Sect. [5] 
we round off with conclusions in Sect. [6] 

Tables [T]and [2] summarise the imaging and spectroscopic data 
presented in this work. A detailed description of the obser- 
vations and data reduction is given in Appendix [A] Tables [3] 
and [4] show the derived aperture photometry. The slit orien- 
tation of the spectroscopic data is illustrated in Fig. Q] for the 
optical spectra and in Fig.|3]for the infrared data. 
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Fig. 1. ADONIS //band image showing the location of the five 
components in the HD 155448 system. For the first time the 
B companion is resolved into two sources, Bl and B2. On top 
of this, we indicate the orientation of the EFOSC2 slits. Image 
orientation here and in all following images is north up and east 
left. 



Fig. 2. EFOSC2 R band image. The cut levels were adjusted to 
also display the arc northeast of the C component, which may 
give the wrong impression that components A and B were not 
resolved. See Fig.Q]for an identification of the components. 



2. Stellar properties derived from optical data 

2.1. HD 1 55448 system members 

In Table [3] we give BVR photometry of the individual com- 
ponents, apart from the Bl and B2 binary that was not re- 
solved at optical wavelengths. Although HD 155448 was of- 
ten referred to as a quadruple system, there is no clear proof 
of the membership of all components to this system. Could 
there be even more components among the neighbouring stars 
in the EFOSC2 field? Therefore, we plot all objects in a colour- 
magnitude diagram (V vs. B - V; see Fig.|Ui. We lack U band 
data for an explicit extinction correction, so we did not dered- 
den the data, but the resulting diagram is clearly different from 
an HRD, thus demonstrating that these stars are not at simi- 
lar distances. A group of seven objects is separated from the 
broader distribution of field stars and appears to be at a simi- 
lar distance. Four of these stars are HD 155448 A, B, C, and 
D. To account for the additional reddening of HD 155448 C 
caused by the CS material, we dereddened this component us- 
ing Ay = 3.7, a value chosen such that the remaining Ay = 
1.5 matches the mean reddening of the other components A, 
Bl, B2 and D (see Table|5]l. We derived the corresponding Ag 
of HD 155448 C according to the interstellar extinction curve 
given in Savage & Mathis ( 119791 1. The remaining three stars of 
the subgroup shown in the upper left of Fig. [4] are located 1— 
2' away from HD 155448. Although they could be at a similar 
distance, they are separated too far to be part of the HD 155448 
system (see also Sect. 12. 51 for a discussion of the distance of 
the HD 155448 system). While we can exclude more system 
members by this test, the membership of HD 155448 C is not 



yet proven, although very likely: it is very unlikely to find five 
B-type stars of luminosity class V on such a small scale (~60 
arcsec 2 ) by coincidence. 

To estimate the likelihood of a chance grouping we consid- 
ered galactic models. However, since the HD 155448 sys- 
tem is located in the direction towards the galactic centre 
(1~353, b — 1-04), we cannot apply common galactic models. 
Alternatively, we investigated the average distribution of BOV - 
B9V stars in this region in the SIMBAD database, and find two 
B-type main sequence stars per deg 2 averaged over a 5 x 5 deg 2 
field, centred on the position of HD 155448. However, this cen- 
sus is incomplete due to the uncertain limiting magnitude in the 
SIMBAD database. More than 90% of the B-type stars found 
in this 5x5 deg 2 field are from the HD catalogue, which ac- 
cording to Flynn & Freeman (1993) is complete down to a 
magnitude V = 9.2 mag (but some objects in our SIMBAD test 
field are as faint as V = 11.5 mag). The number of those B- 
type main sequence stars in the test field brighter than V = 9.2 
mag is one per deg 2 . We want to know how many B stars we 
can expect down to V = 12.2 mag, corresponding to the faintest 
HD 155448 member, if we do not consider component C be- 
cause of its extra extinction. In a standard assumption the num- 
ber of stars multiplies on average by a factor ~2.5 for each 
magnitude (see also, e.g., Table 3 in Roach [19641 1. Based on 
the count of one B-type star brighter than V = 9.2 mag per deg 2 , 
this would result in about 16 B-type stars with luminosity class 
V brighter than V = 12.2 mag that could be expected per deg 2 . 
The HD 155448 components cover an area of ~60 arcsec 2 . 
Therefore, the probability that the entire HD 155448 system 
is a chance alignment is < 10~ 5 . We return to a discussion of 
the membership status in Sect.|4] 
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Fig. 3. VISIR and Spitzer IRS slits and their orientation, plotted 
on top of a SOFI Ks band image. 



2.2. Derivation of spectral types 

To derive the spectral type and luminosity class of the 
HD 155448 candidate members, we carefully compared our 
EFOSC2 spectra with high-resolution (R-80000) spectra 
from the UVES Paranal Observatory Project spectral librarj \ 
(Bagnulo et al. 120031 . the ELODIE (R-42000) spectral librar£ 
(Soubiran et al. 1998), and the theoretical high-resolution 
(R-500000) spectral library BLUEREE0 by Bertone et al. 



(2008). We proceed in a similar way to Carmona et al. (2010). 
First, the template spectra were degraded down to the EFOSC2 
resolution. We then used a dedicated interactive IDL-based 
softwarq^] that allowed us to compare the normalised target 
spectrum simultaneously (i.e., over-plotted) with a normalised 
template spectrum at a user-selected spectral feature. Target 
and template spectra were normalised either by the median or 
by a polynomial fit of second order to the flux in the windows 
used for the spectral comparison. 

We compared our EFOSC2 spectra with the spectral templates 
spectral feature by spectral feature, employing a window of 20- 
100 A width around the central wavelength of each feature. The 
selected spectral features used for comparison are taken from 
the "Atlas of stellar spectra" by Ginestet et al. (1992 ) and the 
features recommended by Morgan, Keenan, & Kellman ( 1943 ) 
from 'An atlas of stellar spectra with an outline of spectral clas- 
sification". To determine the luminosity class, we mainly em- 
ployed the diagnostics of Ginestet et al. (1992 ), described in 
the plates that show the effect of the luminosity for each sub- 
spectral type. In summary, the procedure consisted in the use 
of the Ginestet et al. (1992) and Morgan, Keenan, & Kellman 
( 119431 ) works to guide the choice of which spectral diagnos- 
tics to analyse, and to use the UVES, ELODIE, or BLUERED 
spectral library templates to compare the strengths and shapes 



1 http://www.sc.eso.org/santiago/uvespop/ 

2 http://atlas.obs-hp.fr/elodie/ 

3 http://www.inaoep.mx/~modelos/bluered/bluered.html 

4 Available from A. Carmona: andres.carmona@unige.ch 



8 r 

10 - 
12 - 
14 - 
16 - 
18 - 
20 - 
22 L 



CCDM A 
. CCDM B 



D o 
B 



GSC 




B-V 



Fig. 4. Colour-magnitude diagram of all ~2000 objects in the 
5.5'x5.5' EFOSC2 field. The vertical line in the top left 
quadrant separates stars at the distance of HD 155448 from 
background objects. Letters ABCD denote HD 155448 com- 
ponents (see in the text which assumption was made for 
the C component). CCDM represents the binary star CCDM 
J17130-3216AB. GSC denotes GSC-2 S22231309. Sources 
with B > 21 mag are not displayed. 



of the lines to our spectra. Visually it is relatively straightfor- 
ward to determine which spectral template matches the spec- 
trum under analysis best. However, in cases of degeneracy (es- 
pecially due to rotation), the template spectrum that exhibited 
the smallest residuals (square root of the summed square of the 
difference between the normalised template and target spec- 
trum) was adopted as best match of the target spectrum. 

In Table[5]we present the derived spectral types. We find that all 
objects have spectral type B. In Fig. [5] we show the EFOSC2 
spectra, together with BLUERED synthetic spectra with T e g 
and log(g) of the derived spectral type. In the following para- 
graphs, we describe in detail the spectral classification of each 
source. 



HD 155448 A: The presence of He i in absorption lines and 
the lack of He n indicates that HD 155448 A has a spectral type 
B. The lack of the He n line at 4686 A indicates that the star has 
a spectral type later than BO. Since the Mg n line at 448 1 A is 
much weaker than the He i line at 447 1 A, the spectral type of 
HD 155448 A should be earlier than B5. The presence of the 
Sim lines at 4553, 4568, and 4575 A, the Nn line at 4631 A, 
the O ii lines at 4639 and 4642 A, and the C m lines at 4650 A 
indicate that HD 155448 A has a spectral type earlier than Bl 
if it has luminosity class V, and spectral type earlier than B2 
if it is a giant. The strength of the O n, N n, and C m lines is 
too weak, and the width of the He i lines and hydrogen lines 
(HyS) is too broad to be consistent with luminosity classes I and 
II. While the strength of O n, N n, and C m lines is consistent 
with luminosity classes III and V, the width of the Hei lines 
and H/3 is broader than observed in stars of luminosity class III 
of spectral types Bl to B2. Therefore, HD 155448 A has a lu- 
minosity class V. Altogether, HD 155448 A has a spectral type 
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Table 3. Visual to near-IR photometry of the five components of the HD 155448 system. 



Filter 


A 


Bl 


B2 


C 


D 


Instrument 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 




ID 
L> 


y.l ± U.l 




U. 1 


1 / . J ± U. 1 


1 l.oU ± U.UJ 




V 




IZ.^f ± 


U. 1 


10.1 ± U. 1 


1 l.OU ± U.Uj 


CFVJoV_Z 


K 


( 1 i 
(1) 


IT) , 
1Z.Z + 


U. 1 


1 ^ a _i_ n 1 
1 J.U ± U. 1 


1 1 .4D ± U.UD 


JirUol_.Z 


i 

J 


o. / ± U.l 


11. Q + U.l 


w 


11 a _i_ n i 

ii.y ± u. i 


11 1 _i_ n i 
11.1 ± U.l 


end 
oUrl 


H 


8.8 + 0.1 


11.7 + 0.1 


(2) 


10.8 ±0.1 


11.1 ±0.1 


SOFI 


NB.1.64 


8.5 ±0.1 


11.6 + 0.1 


13.4 + 0.1 


10.7 ±0.1 


11.0 ± 0.1 


NACO 


K 


8.7 + 0.1 


11.7 + 0.1 


(2) 


9.7 ± 0.1 


11.1 ±0.1 


SOFI 


NB.2.12 


8.5 + 0.1 


11.5 + 0.1 


13.1 ±0.1 


9.4 ±0.1 


10.9 ±0.1 


NACO 



Components Bl and B2 cannot be resolved in the visual. Footnotes: (1) saturated, (2) not possible to centre a photometric aperture on B2 due 
to spatial resolution. Please see Appendix lAlf or details about the photometry and its restrictions. 



B1V. We checked several values of t/sin(z), while the model 
spectrum giving the best fit suggested a u sin(/) = 90+5 km/s. 
The error represents the range of v sin(/) values that provide a 
comparable fit to the line profile shapes. The mean value given 
is the one displaying the smallest residuals, but there are a few 
values that are consistent with the observed spectrum. 



HD 155448 B (B1+B2): Because the components Bl and B2 
are spatially unresolved, we consider them first as a single 
source. Given that component Bl is at least one magnitude 
brighter than component B2, the spectrum of HD 155448 B 
is dominated by Bl. Since Hei lines are present at 4471 and 
4923 A, and no Hen lines are observed, HD 155448 B should 
have a spectral type B. Very few Hei lines in the range 4450- 
5100 A are present. Therefore, HD 155448 B should be a late 
B-type star. The Hei line at 4471 A is just slightly stronger 
than the Mg n line at 4481 A. This indicates that HD 155448 B 
should have an earlier spectral type than B7 but later than B4. 
The (weak) strength of the He i lines at 4923 A and the barely 
visible He i line at 4713 and 5015 A indicate that HD 155448 B 
has spectral type B6 or later. Therefore, HD 155448 B most 
likely has the spectral type B6. The broad wings of the HyS 
line are not consistent with the luminosity classes I, II, III, and 
IV, even when considering a vs\n(f) as high as 300 km/s. In 
consequence, HD 155448 B is a dwarf star. We conclude that 
HD 155448 B is a star of spectral type B6V. The fl~3000 model 
spectrum with the best fit suggests a u sin(/) = 150±50 km/s. 

The light of HD 155448 B is dominated by Bl. The near- 
IR images show that the B2 component is 1-2 magnitudes 
fainter than Bl. In the EFOSC2 spectrum, we notice that the 
wings of the H/3 line and the Hei line at 4471 A are slightly 
broader than expected for a B6V star even when including a 
large i/sin(/) (~200 km/s). These spectral characteristics allow 
us to set constraints on the spectral type of the B2 compo- 
nent. Assuming the same distance towards Bl and B2 - and 
the same extinction - the difference in brightness of 1-2 mag- 
nitudes indicates that B2 can have a spectral type ranging from 



Table 4. Mid-IR photometry of HD 155448 C. 



Filter 


C 


C + Arc 


Instrument 




[Jy] 


[Jy] 




PAH2 (11.25 ^m) 


0.51 ± 10% 


0.92 ± 10% 


VISIR 


PAH2_ref (11.88 yum) 


0.44 ± 10% 


0.55 ± 10% 


VISIR 


Q2 (18.72 yum) 


0.73 ± 10% 


1.46 ± 10% 


VISIR 



In the mid-IR narrow-band filters, only the component C is seen. 
Please see Appendix [A3T] for details about the photometry and its 
restrictions. 



B8 to A5. Employing the BLUERED templates degraded at the 
EFOSC2 resolution, we constructed synthetic binary star spec- 
tra with a B6V primary and a B8V to A5V secondary, taking 
care of scaling the fluxes of each component conforming to 
their absolute magnitude. Using our software for spectral com- 
parison, we visually compared the synthetic binary spectra and 
our EFOSC2 spectrum of HD 155448 B. We found that the 
spectrum B6V+B8V is ruled out because the resulting Mgn 
line at 4481 A is stronger than the Hei line at 4471 A. The 
spectra B6V+A1V to B6V+A5V were ruled out because the 
secondary is too faint to have a noticeable effect on the wings 
of the H/3 line. We found that the combinations B6V+B9V 
and B6V+A0V offered the best fits to the H/3 line and the 
relative strengths of the Mgn and Hei lines. As the smallest 
residuals were displayed by the B6V+B9V binary spectrum, 
we suggest that the spectral type of the B2 component is B9V. 
Several values of usin(z) between 100 and 200 km/s gave sim- 
ilar solutions, so we assumed for component B2 a v sin(/) of 
150±50 km/s as derived for component B 1 . 

HD 155448 C exhibits Hei lines in absorption but no Hen 
lines. Therefore, HD 155448 C has a spectral type B. As the 
Hei lines have an FWHM of -240 km/s, HD 155488 C is 
likely rotating fast, at least at ysin(/) of 150 km/s. Emission 
lines are observed most notably Ha, a narrow emission com- 
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Fig. 5. EFOSC2 spectra of the HD 155448 components (black) and rotationally broadened BLUERED synthetic spectra (orange) 
with r e ff and log(g) of the spectral type found. The spectral resolution of the observed spectrum and its model is indicated in 
each panel. The absorption features observed from 6250 to 6350 A are telluric absorption lines. 



ponent inside the H/3 line, and [Oi], [Nn], and [S n] forbidden 
emission (see Table [6}. The lack of a strong He i line in ab- 
sorption at 4713 A indicates that HD 155448 C should be later 
than Bl. At our S IN there is no evidence of strong Sim lines 
at 4553, 4568, and 4575 A, and C in lines at 4647 and 4651 A. 
Therefore, HD 155448 C should have a spectral type later than 
B2. As the Mg n at 448 1 A line is weaker than the He i line at 
4471 A, HD 155448 C has a spectral type earlier than B7. In 
addition, the relative strength between these two lines indicates 
a spectral type later than B3. The strength of the He i, Mg n, and 
the Hy6 line are best matched with the spectral type B4. The lack 
of Si in lines at 4553 and 4568 A and the broad Hfi line wings 
rule out the luminosity classes I, II, and III. The lack of a strong 
He i line at 5047 A is not consistent with luminosity classes III 
and IV. Considering all these characteristics, we conclude that 
HD 155448 C has a spectral type B4Ve. The model exhibiting 
the best fit suggests a v sin(/) = 50±10 km/s. 



HD 155448 D: The spectra of HD 155448 D display broad 
He i lines (v sin(;') ~ 370 km/s). Given that no He n absorption 



lines are present, HD 155448 D has a spectral type B. The He i 
line at 447 1 A and the Mg n line at 448 1 A have a very similar 
strength (with Mgn just slightly stronger). This indicates that 
HD 155448 D has a spectral type later than B6 but earlier than 
B8. The strength of the Hei lines at 4923 and 5015 A is too 
weak to be consistent with the spectral types B6 and B7, thus 
HD 155448 D should have a spectral type B8. Taking the ro- 
tational broadening into account, the width of the wings of the 
H/3 line is too large to be consistent with luminosity classes I, 
II, EH, and IV. As the luminosity class V is the one that matches 
the observed spectrum best, we conclude on a spectral type 
B8V for HD 155448 D. The fl~3000 model with the best fit 
suggests a vsin(/) = 300±15 km/s. 



2.3. Emission lines in HD 155448 C 

Emission lines are only seen in the C component. Their in- 
tensity, in general, is stronger for the slit orientation cover- 
ing the full extent of the arc, slit B-C (see Fig.[TJ, indicating 
that the circumstellar matter to the northeast is also a source of 
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Table 5. Derived spectral types of the HD 155448 compo- 
nents, their corresponding absolute magnitudes, colours, and 
estimated distance ranges. See Sect. 12. 51 for details on the de- 
termination of A v and d. 



Star Type M v (B - V) A v d (kpc) 



A fl 


B1V 


-3.2±0.6 


-0.27= 


t0.02 


1.0+0.5 


1 7 +0 - 7 

1 - / -0.5 


Bl* 


B6V 


-0.9±0.3 


-0.15= 


t0.02 


2.3+0.9 




B2 C 


B9V 


+0.2+0.5 


-0.08= 


t0.05 


2.3+0.9 


1 9 +12 

l y -0.7 


C 


B4Ve 


-1.4+0.2 


-0.19= 


h0.02 


5.2+0.5 


T S +0.7 


D 


B8V 


-0.3+0.4 


-0.11= 


h0.02 


0.6+0.5 


i o +0.6 
l.»_0. 4 



Notes: The errors given for M v and (B - V)o correspond to an error 
in the spectral type of one subclass. 

" Given that the V band image of HD 155448 A is saturated, to 
estimate V, we proceeded as described in Sect. 12.51 
*■ c Since the components B 1 and B2 are unresolved in the B and V 
images, we derived B and V from the observed combined magnitude 
as described in Sect. 12.51 



the emission lines. Ha (single-peaked) represents the strongest 
line, while the forbidden lines [S n], [N n], and [O i] are fainter. 
In the following analysis of the emission lines and their spatial 
extension, we focus on the highest resolution EFOSC2 data, 
observed in August 2008 with Grism 20 in an east-west slit 
orientation, because a contamination by the brighter compo- 
nent cannot be excluded for the A-C slit orientation. 

As v ra d is unknown, the spectrum is centred such that Hei 
at 6678 A - the only photospheric line in Grism 20 - lies at 
v = km/s. This results in a velocity shift of (v ra d + v©) = 
+28 km/s. The error in the fit of a Gaussian to the He i line 
was 9 km/s. Together with the general wavelength calibration 
uncertainty of +4 km/s, the absolute error in the line position 
is then V4 2 + 9 2 ~ ±10 km/s, while the relative error between 
two lines is -y/4 2 + x 2 , with a typical error x m the line centre 
of <2 km/s, resulting in a relative error of ~4 km/s. 

To analyse the spatial extension of the emission lines in the 
eastern direction, we calculated position-velocity diagrams, 
from which we subtracted the continuum PSF (see Fig.|7]i. The 
derived line velocities, widths, and their spatial offset and ex- 
tension in the eastern direction are given in Table [6] We sum- 
marise our analysis of the HD 155448 C emission lines as fol- 
lows. 

- All emission lines are spatially extended in the eastern di- 
rection up to ~4". The [Oi] line at 6300 A is significantly 
less extended than the other lines, while the [Oi] line at 
6363 A is unresolved. 

- All emission lines appear to have a similar velocity shift of 
~10 km/s, except [O i] at 6300 A, which is around ~0 km/s. 
Given that the error on the relative position of the emis- 
sion lines is less than the observed shift (see details above), 
this shift in the forbidden emission lines is most likely real. 
Within their errors, [S n], [N n], and Ha are at the same ve- 
locity. 
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Fig. 6. PMS and main sequence evolutionary tracks according 
to Siess et al. (2000) for initial stellar masses corresponding to 
the spectral types of HD 155448 Bl, B2, C, and D. Overplotted 
are My and colours of all components (red crosses, including 
error bars). The dashed line represents the zero age main se- 
quence. See Sect. l2.4l for further explanations. 



- The [O i] line at 6363 A shows a large error in its velocity, 
caused by the weakness of this line. We do not discuss this 
line in the following, and often refer to the [O i] 6300 A line 
as [Oi] instead. 

The analysis of the line velocities and spatial offsets suggests 
two components in the circumstellar material. One extended 
component, traced by the [Nn], [Sn], and Ha lines, may be 
receding from us, possibly in the form of a wind or outflow 
towards the arc, but a full interpretation of the system geometry 
is difficult (see the end of Sect. l5.3b . A second, less extended 
component, traced by the [O i] line, is close to km/s, but a 
motion in our direction may also be possible within the error 
(cf. Table [6). As the projected velocities are much lower than 
100 km/s, which would be a characteristic value for a massive 
outflow or jet, the gas may be moving in a direction close to the 
plane of the sky, but the low velocities could also be intrinsic. 
We return to a discussion of the forbidden emission lines in 
Sect.E2 

Could the arc matter have influenced the spectrum of 
HD 155448 C? As shown in Fig.|7]and in Tabled we see that 
the peaks of most emission lines are inside a region correspond- 
ing to the PSF of HD 155448 C. At the given spatial resolution, 
this indicates that most of this emission is associated with the 
star. If the forbidden emission was dominated by the arc mat- 
ter, we would have observed a clearly higher spatial offset of 
the emission line peaks. Nevertheless, some contribution of the 
arc matter can be expected as well. In particular, the arc's inner 
border towards HD 155448 C cannot be defined clearly, and 
the forbidden emission lines seem to originate in a transition 
region. 

High spatial and spectral resolution spectroscopy is required to 
determine the contribution of arc matter and to better constrain 
the velocities of the gas. 
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Table 6. Emission lines observed in HD 155448 C for the east-west slit orientation. 



1 

A 


Line 




r W rllVl 


Centre 


vjnsei 


Extension in E 


t A\ 

(A) 




(A) 


LKm/sj 


[Km/sj 


L J 


L J 


6300 


[Oi] 


-1.0 


100 ±6 


-4+ 10 


-0.24 ± 0.06 


-2.0 


6363 


[Oi] 


-0.3 


100 ± 19 


+ 13 ± 14 


-0.24 ± 0.06 


unresolved 


6548 


[Nil] 


-1.7 


91 ±6 


+ 13 ± 10 


-0.60 ± 0.06 


-3.0 


6563 


Ho- 


-23 


95 ±4 


+ 13 ± 10 


-0.60 ± 0.06 


-4.5 


6583 


[Nil] 


-4.3 


91 ±4 


+ 11 ± 10 


-0.84 ± 0.06 


-3.5 


6716 


[Sn] 


-4.0 


88 ±4 


+9± 10 


-0.36 ± 0.06 


-3.5 


6731 


[Sn] 


-6.1 


88 ±4 


+ 10+ 10 


-0.36 ± 0.06 


-3.5 



The error of the FWHM is given by the relative wavelength calibration uncertainty and the uncertainty of the Gaussian fit. For the line centres, 
the error includes the absolute wavelength calibration uncertainty (see Sect. l2.3t . The spatial extension was determined via position-velocity 
diagrams (see Sect. |2.3l >. 



2.4. Age 

In Fig. [6] we plot the derived My and colours from Table [5] 
in a colour-magnitude diagram with PMS and main sequence 
evolution tracks from Siess et al. (2000). Only the tracks 
for the initial masses corresponding to spectral types of the 
HD 155448 system are shown, which is 3.0 M© (B2), 3.5 Mq 
(D), 5 Mq (Bl), and 6 Mq (C). Because the Siess models are 
only available up to 7 Mq, corresponding to spectral type B3, 
HD 155448 A would be located above these tracks, as its mass 
is higher than 7 Mq. Within the errors, all stars are very close 
to the ZAMS, thus showing that the system is not a post-AGB 
or pre-PN object, as stated in the literature. 



2.5. Distance 

The error on the Hipparcos parallax (1.65+1.99 mas; Perryman 
et al. 1 19971 1 is too large to derive a meaningful distance, which 
nominally would correspond to 606 *™j pc, with an open up- 
per distance limit. Therefore, to constrain the distance of the 
HD 155448 components, we have to rely on the derived spec- 
tral types and photometry. In Table|5] we list the appropriate ab- 
solute magnitudes for the corresponding spectral classification, 
obtained from Schmidt-Kaler (|19821 l. The estimated reddening 
Ay was derived from the observed B - V colours, the intrin- 
sic colours given by Schmidt-Kaler (1982), and it assumes a 
standard reddening law (R = 3.14) 

A v - 3.14 [(B -V) - (B - V)o] ■ 

The distance follows from the classical distance modulus: 
m v - My = 5 log(fi) - 5 +A V . 

Since the V band image is saturated for component 
HD 155448 A, to estimate V for component A we used B - V = 
0.06 as given in SIMBAD, and B = 9.05 derived from our 
EFOSC2 data. Using the Hipparcos B —V = 0.186 would re- 
sult in a ~0.3 kpc closer distance for HD 155448 A. The com- 
ponents Bl and B2 are not resolved in the EFOSC2 images, 
but only separated in the infrared. We estimate their individ- 



ual B and V magnitudes from the observed combined mag- 
nitude, by using the difference in the expected absolute mag- 
nitudes between spectral types B6V and B9V (AMy = 1.10, 
AMb = 1.18). For example, to derive the V magnitude of each 
component we used 

V B i = Vbi +b2 + 2.5 log [1 + ioCA*v/-zS] 

V B 2 = VB1+B2 + 2.5 log [1 + 10 C^v/2.5)] _ 

In this way, we obtained V = 12.7 ± 0.2 mag and B = 13.3 ± 
0.2 mag for the component Bl, and for component B2, V = 
13.8 ± 0.2 mag and B = 14.5 ± 0.2 mag. 

We summarise the results of the distance determination in 
Table [5] The distances to all components are consistent within 
their errors, with a mean distance of 2.2 kpc. While the uncer- 
tainty in the distance estimate seems large, we emphasise that 
these error ranges are mainly caused by an uncertainty of only 
one subclass in the spectral type derivation. 



3. Results derived from infrared data 

3.1. A fifth component in the HD 155448 system 

To our knowledge, this paper offers the first high-resolution, 
adaptive optics observations of the HD 155448 system. The 
ADONIS data reveal that the component HD 155448 B can 
be resolved into two stars (cf. Fig.[T|i with a separation of 
1.21+0.04". This additional companion is also visible in 
NACO and SOFI data. In NACO data, the separation is 
1.16+0.01", and 1.14+0.05" for SOFI. Subpixel accuracy was 
applied in calculating the component separations. We do not 
want to speculate about whether we see any motion during 
the five years between these data sets, because all results 
(marginally) overlap within their error ranges. The photome- 
try is given in Table [3] 
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Fig. 7. Important emission lines observed in HD 155448 C, extracted from spectra obtained in an east-west slit orientation. The 
upper panel displays the extracted spectral line, the middle panel shows the position-velocity diagrams in the raw 2D spectrum, 
while the lower panels expose the position-velocity diagrams after the subtraction of the continuum PSF. The first (outer) contour 
represents the 5 <x level with respect to the background noise. Subsequent contours are at 10 cr intervals, except for Ha, which is 
composed of 30 cr intervals. 



3.2. Imaging and photometry 

In the images we found arc-shaped emission, curving from 
HD 155448 C towards the northeast. This feature is seen both 
in the near- and mid-IR, and is also visible in optical data. The 
extension is up to ~3" from HD 155448 C in all the observed 
passbands. A difference in length of 0.2" between the filters 
may depend on the sensitivity of the individual data sets. The 
inner end of the arc, facing towards HD 155448 C, is not clearly 
constrained. At optical wavelengths, in the near-IR and in the 
mid-IR around 10 fim, little extra emission is seen very close 
to the star, while significant emission starts at distances beyond 
-1.0" from HD 155448 C. In the mid-IR PAH band filter and 
at shorter wavelengths, the intensity is maximal in the arc's 
centre. The appearance in the Q band filter is clearly different, 
where the arc extends all the way towards the star. Moreover, 
the peak of the emission is now located close to the star, in an 
area that shows nearly no emission in the near-IR or at visual 
wavelengths. 

Tables [3] and H list the photometry for all passbands used in 
this work. Only for adaptive optics data was it possible to ob- 



tain resolved photometry for the B2 companion. The SOFI 
and NACO magnitudes agree well within their errors, despite 
the problems calibrating the NACO data (see Appendix IA.41 i. 
For VISIR we performed PAH-on (11.25 ^m) and PAH-off 
differential imaging (continuum around 11.88 /mi), and sub- 
tracted both frames after a re-centring with 0. 1 pixel accuracy. 
For a description of the errors in the photometry, we refer to 
Appendix fA.5. II In Table 0] we also give for VISIR the total 
flux of HD 155448 C plus the circumstellar matter in the arc, 
because resolving both is not accurately possible in the Q band, 
as is recognisable in Fig. [8] 

3.3. Mid-IR spectroscopy 
3.3.1. VISIR 

We used VISIR low-resolution spectra to study the composition 
of the circumstellar matter close to HD 155448 C and in the arc. 
None of the other stellar components is detected in the mid-IR, 
what was also seen in our previous test observations using the 
TIMMI2 mid-IR camera at ESO La Silla Observatory, which 
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PAH on PAH off PAH only Q2 




Fig. 8. VISIR images in the PAH2 filter at 1 1 .25 yum ("PAH on"), the continuum at 1 1 .88 /xm ("PAH off"), the difference of both 
previous images ("PAH only"), and in the Q2 filter at 18.72 //m. Only the C component is seen in the mid-IR narrow-band filters 
used. The VISIR slit orientation at 50° position angle is indicated in the left image (cf. also Fig. [3]). To obtain optimal contour 
plots for the different intensity ranges, 50 equidistant contour levels were plotted for "PAH on" and "PAH off", 5 contour levels 
for "PAH only", and 20 levels for the Q2 image. 



has a larger field of view than VISIR. However, our mid-IR 
narrow-band filters used cover only a small - although repre- 
sentative - range of the N band. An eventual mid-IR emission 
originating in the remaining HD 155448 components at other 
N band wavelengths cannot be fully excluded. The slit orienta- 
tion and the procedure for disentangling the emission near the 
star from that of the arc is explained in Appendix lA.5.21 



HD 155448 C : The spectrum of HD 155448 C is dominated 
by silicate emission between ~8— 12 /urn (Fig. [9] left panel). 
For a description of silicate emission features at these wave- 
lengths, we refer, e.g., to Schiitz et al. (2009) and references 
therein. The silicate feature of HD 155448 C can be interpreted 
as emission mainly from small, amorphous silicate grains lo- 
cated in the hot surface layer of a circumstellar disk. A smaller 
peak at 11.3 z/m indicates either the presence of crystalline sil- 
icates or of polycyclic aromatic hydrocarbons (PAHs), which 
both exhibit an emission peak around this wavelength. We pro- 
vide more details on the dust composition in Sect. 13.3.21 



Arc : Unlike for the spectrum of the stellar component, we 
see no silicate features related to the arc. Here, the spectrum is 
dominated by strong PAH bands at 8.6, 11.3, and 12.7 yum. We 
also see fainter PAH features at 1 1 .0 and 12.0 fim. In Table|7]we 
give a detailed list of the measured properties of the observed 
features and of the bending modes that can cause them (from 
Hony et al. 1200 11 1. Based on the spectral appearance of PAH 
bands, Peeters et al. (120021 ) define characteristic PAH classes 
A/B/C, which represent different stages of PAH processing. 
By comparing the peak position of the PAH bands in our spec- 
trum to these classes, it is clear that the arc's PAH features do 
not belong to the highest processed class C, as there is a strong 
feature at 8.6 fim that is absent in class C sources (Peeters et al. 
2002). To decide on class A or B, we need to also consider the 
PAH bands between 6-8 yum, and refer to the Spitzer spectrum 
in the next section. In Fig.[l0] we compare the PAH spectrum 



of the arc with that of the mean ISM and that of HD 169142, 
a prototype of PAH features observed in Herbig Ae/Be stars 
(see Acke et al. I2010t . It is clear that the arc spectrum is more 
like that of the ISM than of the Herbig Ae star, especially when 
focussing on the region near 1 1 fim. The small feature at the 
base of the 11.3 fim band towards shorter wavelengths is the 
11.0 fim PAH band, which is attributed to the out-of -plane 
bending of CH in PAH cations. Therefore, its emission is an 
excellent tracer of ionised interstellar PAHs. Indeed, in a large 
sample of isolated Herbig Ae/Be stars, the 11.0 fim feature is 
absent (Acke et al. 1201 Oi l. This suggests that the arc-shaped, 
PAH emitting region is not directly related to the disk around 
HD 155448 C, but rather an ISM remnant of the stellar forma- 
tion process. 



3.3.2. Spitzer 

As described in Sect. 13. 3. Tl we see the emission from 
HD 155448 C and the arc in the mid-IR. However, in the 
Spitzer spectra the emission from both sources is not spatially 
resolved. To derive the composition and physical parameters 
of the dust grains, whose emission features are clearly visible 
at mid-IR wavelengths, we used the two-layer temperature dis- 
tribution (TLTD) method. For details of the spectral analysis 
method, as well as for the applied optical constants of the dust 
species used, we refer to Juhasz et al. (2009). The spectrum was 
fitted in two wavelength regions: short (5-17 /um) and long(17- 
35 fim), with five dust components: amorphous silicates with 
olivine and pyroxene stoichiometry, crystalline forsterite and 
enstatite, and amorphous silica. Three different grain sizes were 
used for the amorphous silicate and silica components (0.1, 1.5, 
and 6.0 fim), while only two for the crystals (0.1 and 1.5 /um), 
as no evidence was found for large crystals in the spectrum. A 
mean ISM PAH template from Hony et al. (2001) was used to 
derive the individual PAH band profiles for the fits. The results 
are shown in Fig. [TT] and Table [8] 
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HD 155448 C (without the arc) Arc (without HD 155448 C) 
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Fig. 9. VISIR N band spectra. Left panel: HD 155448 C, showing a silicate feature. Right panel: the arc, with strong PAH bands 
at 8.6, 11.3, and 12.7 fim. The y-axis scale is not identical in both figures, but was adjusted to optimally show the features. The 
gap between approximately 9.5 fim and 10.3 fim is due to incomplete instrument coverage (see Appendix lA.5.21 i. 



Silicates : The absolute flux level in the Spitzer IRS spectrum 
is higher in the 8-13 fim interval than in the VISIR spectrum, 
indicating that the Spitzer spectra may include more emis- 
sion from the surrounding material due to the wider slit. The 
dust composition in the 10 fim region is dominated by small 
(0.1 fim) amorphous grains. The low peak-to-continuum ratio 
of the 10 fim silicate feature, which is usually explained by the 
presence of large grains, is probably caused by the high extinc- 
tion towards HD 155448 C. Crystalline silicates represent only 
about 13 % of the total silicate mass in the short interval. This 
is a normal value among Herbig Ae/Be stars (see, e.g., Juhasz 
et al. 120 101 or van Boekel et al. 2005). The crystal population 
is dominated by forsterite, whereas the contribution of enstatite 
is negligible. When comparing the fit of the short interval with 
that of the long interval, we see that the forsterite-to-enstatite 
mass ratio decreases with wavelength (i.e. with radius in the 
disk), similar to what is observed in other Herbig Ae/Be sys- 
tems (Juhasz et al. 120101 ), and even T Tauri stars (Meeus et 
al. 12009). We note that, although the general quality of the fit 
is good, there are mismatches between the observed and the 
model spectrum. In the short region these mismatches are re- 
lated to the PAH features, because their shape was not perfectly 
reproduced by our model templates. In the long wavelength in- 
terval, the crystalline silicate bands in our model seem to be 
systematically shifted towards longer wavelengths than in the 
observed positions. These uncertainties in band positions are 
likely related to grain shape effects and/or to the dependence 
of the mass-absorption coefficients on the temperature that is 
not included in our analysis (see also Juhasz et al. 120101 for a 
discussion on this topic). 



PAH features : The observed shape and positions of the PAH 
bands in the spectrum of HD 155448 C are consistent with that 
of a Class A profile in Peeters et al. ([2002), typically found in 
the ISM (cf. also Sect. l3.3.fl paragraph on the arc). 



- The dominance of the 7.6 fim sub-band over that of the 
7.8 fim sub-band in the 7.7 fim PAH feature complex sug- 
gests that the emitting molecules are small (number of car- 
bon atoms Nc < 90-100) and ionised (Bauschlicher et al. 
120081 . 

- Acke et al. {2010) find decreasing aliphatic over aromatic 
line strength ratios with increasing stellar effective temper- 
ature. This correlation is explained by the higher stability 
of the aromatic bonds against the destructive UV photons 
compared to the aliphatic bonds. Given that no evidence 
was found for aliphatic emission features (at 6.8 fim and 
7.2 fim) in its spectrum, HD 155448 seems to follow this 
trend. 

- Hony et al. ( 120011) and Keller et al. ( 120081 ) interpret the 
6.2/1 1.2 fim and 12.7/1 1.2 fim line strength ratios as indica- 
tors of the ionisation state and irregularity of the molecules. 
Comparing the calculated line strength ratios of HD 155448 
(F6.2/F11.2 = 1.52 + 0.12, F12.7/F11.2 = 0.42 + 0.04) to 
Fig. 5 in Hony et al. (2001), the PAH molecules in the 
arc near HD 155448 C are probably ionised and irregular, 
pointing to PAH processing by the relative intensive UV 
field in the HD 155448 environment. 



4. Position and motion 

We do not aim to give absolute astrometry, since the coordi- 
nates were not derived by means of astrometric standard stars. 
At the time of our observations, both ADONIS and EFOSC2 
data had not yet been offered with a world coordinate system, 
which we were able to fix for the EFOSC2 files. Obviously, the 
NACO adaptive optics data permit a much more precise cal- 
culation of positions and angles than the SOFI frames, but the 
rather short time in between these two data sets (one year) does 
not permit calculating component motions at the distance of the 
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Table 7. Properties of the PAH features. 



Feature 


Bending Mode 


Peak Position 


Central Wavelength 


(jim) 




(/an) 


(pm) 


6.2 


CC stretching 


6.21 


6.23 


7.6 


CC stretching 


7.58 


7.58 


7.8 


CC stretching 


7.82 


7.85 


8.6 


CH in plane 


8.64 


8.66 


11.0 


CH out of plane 


11.03 


11.02 


11.3 


CH out of plane 


11.25 


11.26 


12.0 


CH out of plane 


12.01 


12.01 


12.7 


CH out of plane 


12.75 


12.63 



Table 8. Dust composition of the Spitzer IRS spectrum, in two 
wavelength regions. 



The central wavelength was determined by calculating the weighted 
centre of the feature. 




Fig. 10. Comparison of the PAH bands in the mean ISM (from 
Hony et al. l250TT l. in HD 169142 (IRS spectrum), and in the arc 
near HD 155448 C. The feature at 1 1 .0 //m clearly differentiates 
the arc matter near HD 155448 C from HD 169142, while the 
feature is also present in the ISM. 



HD 155448 system. Therefore, we focus only on the NACO 
data in Table [9] The errors in our RA and Dec can be deter- 
mined from the positional uncertainty in the data sets of various 
nights, which were taken with different guide stars. (The astro- 
metric accuracy is limited by the accuracy of the guide star 
coordinates.) We do not quote RA and Dec for the literature 
data, as the VizieR web interface points out that their accuracy 
could be low. 

Can we see evidence for a gravitationally bound motion 
of the HD 155448 components? The proper motion of 
(5RA = 2.77 mas/yr and £Dec = -0.85 mas/yr, together with 
their uncertainties of the same order (Perryman et al. 119971 the 
catalogue entry refers only to HD 155448 A), could - during 
the past ~80 years - have caused a change in p up to ~0.4" 
for unbound stars, depending on how the mean proper motion 



Parameter 


Grain 


Mass fraction/ 


Mass fraction/ 




size [/mi] 


PAH intensity 
(5-17/zm) 


PAH intensity 
(17-35jum) 


X' 




61.8 


91.1 


Amorph. Olivine 


0.1 


71 2 +2 - 2 


32 1 +3 - 7 


Amorph. Olivine 


1.5 


8.2 + } 2 \ 


53 1 +43 


Amorph. Olivine 


6.0 






Amorph. Pyroxene 


0.1 






Amorph. Pyroxene 


1.5 






Amorph. Pyroxene 


6.0 






Cryst. Forsterite 


0.1 


2.8 + _° j 


c /"» +0 3 

5 2 

-0.3 


fr\7ci~ T-ir\i"c t"F*t*i tc* 
V^IVaL. .TLflalClllC 


1 c; 
i .j 


i n 3 +0 - 6 




Cryst. Enstatite 


0.1 




I O +0.2 
L,y -0.2 


Cryst. Enstatite 


1.5 






Amorph. Silica 


0.1 


5 3 +03 

J - J -0.3 




Amorph. Silica 


1.5 


1 o +0 - 8 

l,u -0.6 


5 3 +0 - 2 

' -1.7 


Amorph. Silica 


6.0 




2 3 +L2 


PAH 6.2 
PAH 7.6 
PAH 7.8 
PAH 8.6 
PAH 11.3 
PAH 12.7 




0.285 ™» 

n 997 +0.005 
-0.004 

189 +0 003 
U.105- -0 Q02 

226 +0 003 

-0.003 

r\ tc9 +0.006 
-0007 

fl 1 ~\1 +0 003 
u - u ' -0.003 





Relative silicate mass fractions and their errors are given in percent. 
The mass fractions marked as non-detections are either fitted as 0% 
or have < l<r significance. The numbers indicated for PAHs are 
relative intensities compared to the normalised PAH model profile 
from Hony et al. OTTO . 



and its uncertainty would coadd. This value could eventually be 
doubled, when the random motion of two stars would be in the 
opposite direction. Such a large motion is not seen between the 
components A, B, and C. On the other hand, the mean proper 
motion and its uncertainties could also result in lower values. 
Without knowledge of the precise coordinates from the years 
1911 and 1928, respectively, we cannot analyse the precise tra- 
jectories of the components. However, no significant angular 
motion is seen between now and ~80 years ago. 



5. Discussion 

5.1. SED of the HD 155448 system 

Our analysis shows that HD 155448 is a young quintuple sys- 
tem, in contrast to earlier descriptions in the literature refer- 
ring to it as a post-AGB star. This earlier classification was, 
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A [am] A H 

Fig. 11. Fit of the Spitzer IRS spectrum of HD 155448 C. The points in the upper panels represent the Spitzer data, while the 
overplotted continuum line is the fitted model. See Table[8]for the individual contributions of the model dust species. The residuals 
of the fit are shown in the lower panels. The quality of the fit is good in general; however, there are mismatches between our 
model and the observed spectra. In the 5-17 fim region most of the mismatches are related to the strong PAH emission features. 
For the 17-35 fim region the mismatches are instead caused by differences in the peak positions of the strong crystalline silicate 
bands between the model and the observed spectrum. 



however, only based on IRAS colours, which does not permit a 
unique classification (cf. Sect.[TJ. 

Based on the interpretation of the SED, Malfait et al. (1998 ) 
derived the presence of a large inner hole in the disk of 
HD 155448 A. However, due to missing spatial resolution in 
the photometry, their SED merged the fluxes of HD 155448 A, 
dominating the visual wavelengths, and the mid-IR emission 
from HD 155448 C. In Fig.Qj] we plot resolved SEDs for 
both components. Our photometry is dereddened and plotted 
in comparison to a Kurucz atmosphere model (Kurucz l"l9941 ). 
The sum of the spatially resolved photometry of HD 155448 A 
and C agrees well with the spatially unresolved data of Malfait 
et al. ( 119981 1. This means that, instead of a single transitional 
disk, the HD 155448 system hosts four stars with very little 
or no circumstellar material and one with a disk. Substantial 
NIR excess emission is seen for HD 155448 C. If this emission 
indeed originates in a circumstellar disk, there is no reason to 
assume an inner hole, but rather a hot inner rim. 

Besides the quintuple, we uncovered an arc-shaped emission 
region northeast of HD 155448 C. We exclude the possibility 
that the arc could be a chance alignment because in the Q band, 
we see that the circumstellar matter is connected with the star 
(Fig. [8}. But what is the origin of this circumstellar matter and 
what causes the observed elongated curved emission? 



5.2. Origin of the forbidden emission lines 

For late-B and A-type Herbig stars, the formation of a circum- 
stellar disk during the stellar formation process is theoretically 
predicted. The observed geometry of the circumstellar matter is 
less clear for early-B stars, where the material may be present 



both in spherical (envelope) and disk-like geometries. The clas- 
sical Be star phenomenon is explained by ejected gaseous mat- 
ter in an equatorial disk around a rapidly rotating star (see, e.g., 
Porter & Rivinius 2003 ). In this model, line emission other than 
from hydrogen is typically much weaker. Neither the strong 
forbidden lines nor the arc-shaped circumstellar matter could 
be explained with this scenario. 

The forbidden [O i] emission line is observed in a wide range of 
astrophysical objects, and seen around both young and evolved 
stars. For Herbig Ae/Be systems, Acke et al. (12005) find a cor- 
relation between the strength of the [Oi] line at 6300 A and 
the far-IR excess and show that stars with a flared disk geome- 
try are more likely to exhibit the [O I] line than self-shadowed 
disks. In these objects, the [Oi] line observed is centred at the 
star's velocity and is double peaked, therefore indicating that 
the emission arises from the circumstellar disk. The appearance 
of [O i] together with [S n] in emission is commonly observed 
for accretion phenomena and outflows both in massive and low- 
mass stars (e.g., Hartmann & Raymond 1 19891 Hartigan et al. 
119951 ), while both lines are quite common in Herbig Ae/Be sys- 
tems. Vieira et al. (120031 find the forbidden [O i] lines (6300 
and 6363 A) and [S n] lines (67 1 6 and 673 1 A) for 46% of their 
131 Herbig Ae/Be candidate stars, and find a higher occurrence 
of those lines among B-type stars and those objects that have 
stronger circumstellar emission. 

We see indications of a wind or outflow from HD 155448 C. 
The velocities of the [S n] and [N n] lines, which are spatially 
extended in direction towards the arc, point to a wind that is 
receding from us (cf. Sect. l2.3b and that probably may interact 
with remnant matter of the star formation process located in 
the arc. Indeed, jets and collimated outflows from early B-type 
stars are observed for objects on the ZAMS (see, e.g., Shepherd 
120051 ). 
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Table 9. Position angles and separations relative to component HD 155448 A. 



Observation 


Source 


C omponent 




n 
r 




Q 


RA 


Dec 


(year) 








["] 




[°] 


[J2000.0] 


[J2000.0] 


1928 


DN2002 


A 
















1928 


DN2002 


B 




2.9 




280 






1928 


DN2002 


C 




4.2 




16 






1911 


DN2002 


D 




9.9 




91 






2005 


NACO 


A 










17:12:58.79 + 0.01 s 


-32:14:33.7 ± 0.1 s 


2005 


NACO 


Bl 


2.84 


±0.01 


279.8 


± 0.3 


17:12:58.57 ± 0.02 s 


-32:14:33.2 ± 0.1 s 


2005 


NACO 


B2 


2.26 


± 0.01 


302.7 


± 0.3 


17:12:58.64 ±0.01 s 


-32:14:32.5 ± 0.1 s 


2005 


NACO 


C 


4.39 


± 0.01 


16.6 


± 0.3 


17: 12:58.89 ± 0.02 s 


-32:14:29.5 ±0.1 s 


2005 


NACO 


D 


9.43 


± 0.01 


90.4 


± 0.3 


17:12:59.54 ± 0.02 s 


-32:14:33.8 ± 0.1 s 



References: DN2002 (Dommanget & Nys 2002 1, NACO (data in this work). The errors of the position angles and separations in the literature 
data are not known. 



Why is the [O i] emission not part of this receding wind? In a 
study of forbidden emission lines in T Tauri and Herbig Ae/Be 
stars, Hamann et al. (1994) found that the different lines are 
formed at distinct densities and temperatures in the circumstel- 
lar environment. While [Oi] needs a relatively high density, 
this is much lower for [N n] and even lower for [S n]. If an out- 
flow exists from HD 155448 C, it implies a radial decrease in 
density with distance from the star. This agrees with our obser- 
vation, where [O i] is found closer to the star than the other two 
forbidden lines (cf. Table |6). These regions overlap to some 
extent, but each one is dominated by a different transition. 

Multicomponent winds were also found in, e.g., DG Tau 
and HL Tau (Solf & Boehm 1 1993ft , where the low-velocity 
[Oi] 6300 A components coincide with the star, while the 
high-velocity components [N n] and [S n] are distributed 
along the jet direction. No high-velocity components exist in 
HD 155448 C, but within the different velocities, this system 
might be a similar scenario. The low wind velocities may be 
intrinsic or caused by a projection effect. Some of the ob- 
served [Oi] emission might also originate in the disk around 
HD 155448 C, but the link to a wind origin relies on the fact 
that the [O i] line is observed slightly extended to the east (see 
Fig.0. 

We note that the interpretation by means of a wind or outflow is 
not a unique solution of the HD 155448 system geometry, and 
discuss an alternative possibility at the end of Sect. \5~3\ 

5.3. The arc 

We could not find a scenario that perfectly explains the emis- 
sion pattern observed around HD 155448 C. The silicates and 
PAHs represent remnant interstellar matter, the first one proba- 
bly located in a disk around HD 155448 C, the latter in the more 
extended, asymmetric arc. There are two possible scenarios: 

- A circumstellar disk or envelope exists around 
HD 155448 C, with the MIR flux dominated by emission 



from hot silicate grains. The arc to the northeast consti- 
tutes remnant matter from the stellar formation process. 
The PAHs in the arc are excited by the irradiation of 
all close B-type stars. The wind or outflow originating 
in HD 155448 C impacts the ISM and thus creates its 
arc-shaped form. 

- In an alternative scenario we may speculate on the forma- 
tion of another, so far undetected, late-type protostar in the 
HD 155448 system, located at the central part of the arc. As 
the stellar formation timescales are much larger for lower 
masses, typically 10 7 vs. 10 4 years, this protostar could still 
be embedded in its parental cloud. Indeed, in the mid-IR we 
see a knot in the central part of the arc, which is less promi- 
nent in the near-IR and not observed at optical wavelengths. 
The arc's shape may have been caused by the interaction of 
the wind with the envelope of this hypothetical low-mass 
object. 

We favour the first scenario, because, although the other alter- 
native is plausible, we were not able to provide more support 
for it with this work. Altogether, we suggest the following in- 
terpretation: there is an outflow from HD 155448 C in the di- 
rection of the arc. The latter is composed of remnant material 
of the star-forming process and illuminated by the radiation of 
the stars. Its PAHs are partly ionised and resemble ISM ma- 
terial. The forbidden emission lines trace several wind compo- 
nents along the outflow, with [O i] centred on the star, where the 
outflow density is highest and the velocity component close to 
zero, while [N n] and [S n] are formed at greater distances from 
HD 155448 C, where the outflow density is lower and the wind 
velocity higher. 

An important limitation of the outflow interpretation is that 
only the receding (redshifted) outflow component is seen. 
Because in our observations we do not have evidence of ma- 
terial blocking the preceding (blueshifted) outflow component, 
we must conclude that the outflow in HD 155448 C has only 
the single component receding from us. There may be the fol- 
lowing explanations. 
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Fig. 12. Spectral energy distributions of the HD 155448 A and 
C components constructed from spatially resolved photome- 
try. Overplotted are the fitted Kurucz atmosphere models. See 
Sect. l5.1l for an interpretation. 



- Lovelace et al. ( 1201 Oi l show that an intrinsically one-sided 
wind can indeed occur for a rapidly rotating protostar with 
a magnetosphere in the propeller regime and with a mag- 
netic field consisting both of a dipole and a quadrupole 
component. Whether such a configuration may apply to 
HD 155448 C is not known, since additional data would be 
needed for a detailed investigation of the magnetic field. 

- Alternatively, we may also interpret the redshifted forbid- 
den line emission as originating in infalling remnant matter 
of the star-forming process. In this case the arc material will 
be situated in the foreground of HD 155448 C. 

With the data presented in this work, we cannot decide which 
of these explanations is the more likely one. Spectra with both 
higher spectral and spatial resolution, ideally obtained with an 
integral field unit, are needed to better understand the physical 
interactions between HD 155448 C and its surrounding matter. 
To set constraints on a magnetic field of HD 155448 C, spec- 
tropolarimetric observations at a high spatial resolution are re- 
quired. Given the rather faint magnitude of HD 155448 C, these 
follow-up observations will be ambitious. 



6. Conclusions 

We presented multi-wavelength photometry of the HD 155448 
system, covering the optical and near- and mid-IR ranges, from 
the B to the Q bands. From the optical spectra we derived the 
spectral types of all five components. HD 155448 C shows 
strong emission lines, while circumstellar matter in an unusual, 
arc-shaped form is seen within 3" northeast of HD 155448 C. 
The main results of this study follow. 

1 . Near-IR images obtained with ADONIS, SOFI, and NACO 
revealed that HD 155448 B is a binary. Including this dis- 
covery, the system now comprises five objects within ~12". 
Our distance estimate from spectral parallaxes indicates 



that all stars are at a similar distance of ~2 kpc, but the error 
of 0.7 kpc is rather large, stemming from the fact that B star 
luminosities vary strongly with one subclass uncertainty in 
the spectral classification. 

2. From EFOSC2 optical spectroscopy of all five sources, we 
determined their spectral types as B1V (HD 155448 A), 
B6V (HD 155448 Bl), B9V (HD 155448 B2), B4Ve (C), 
and B8V (D). 

3. In a colour-magnitude diagram with stellar evolutionary 
tracks from Siess et al. (12000b . the HD 155448 components 
B 1 , B2, C, and D are located very close to the ZAMS (while 
the mass of HD 155448 A has no corresponding track in 
the Siess models). Earlier classifications of HD 155448 as 
a transition object from post-AGB to pre-PN are not con- 
sistent with this location. 

4. HD 155448 C exhibits strong emission lines of Ha and 
[S n], as well as lines of [Nn] and [Oi]. All emission lines 
in HD 155448 C are spatially extended in the eastern direc- 
tion up to ~4"and appear to have a similar velocity shift, 
except the [O i] line, which is significantly less extended in 
spatial direction and close to the rest velocity of the star. 

5. No significant orbital motion in the HD 155448 system was 
observed, neither in the time span of our data sets between 
the year 2000 to 2005 nor when considering literature data 
from ~80 years ago, though the uncertainties in the latter 
data did not permit an exact comparison. 

6. Only HD 155448 C is detected in the mid-IR. From VISIR 
10 jum spectroscopy, we see that the C component is sur- 
rounded by small, warm silicate grains, suggesting the pres- 
ence of a circumstellar disk. In contrast, the emission fea- 
tures in the arc are entirely of PAHs. Their PAH band sig- 
nature is more characteristic of ionised matter than of a cir- 
cumstellar disk. 

7. We find no scenario that could satisfactorily explain the 
shape and nature of the arc. A wind or outflow originating 
in HD 155448 C may impact the ISM, thereby creating its 
arc-shaped form. Alternatively, we may speculate about an 
interaction of the wind with the envelope of a hypothetical 
low-mass object. 

8. High spatial and spectral resolution observations of gas 
lines are needed to further study and understand this enig- 
matic system. 
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Appendix A: Details on the observations and data 
reduction 

A.1. ADONIS NIR imaging 

The ESO Adaptive Optics Near Infrared System (ADONIS; 
Beuzit et al. 1 1997ab was available on the 3.6 m telescope at 
ESO La Silla until September 2002. Typical Strehl ratios were 
around 0.1 in J band and 0.3 in H band. ADONIS was coupled 
with the near-IR camera SHARPII+ which operates in the / to 
K band. Our data were acquired on June 8, 2000 in the H and 
S K (short K) filter with total exposure times up to 600 seconds. 
We attached a fully opaque coronographic mask in front of the 
Lyot (pupil) stop to reject the peak of the PSF for the brightest 
source, to increase the integration time and sensitivity in order 
to reveal fainter structures. Details on the coronograph and its 
performance are given in Beuzit et al. ( 1997b). For HD 155448, 
a mask of 1 .4" and a lens scale with 0. 1 "/pix were used. We co- 
added ten exposures of 60 seconds in both filters. The mask's 
fixed position meant that no jitter offsets were possible. Data 
cube clean-up, dark and flatfield correction, as well as sky sub- 
traction were performed in standard ways. Figure Q] shows an 
ADONIS H band exposure. No PSF subtraction was applied in 
this case. 

Despite the high sensitivity of the ADONIS data, no reliable 
photometry could be obtained. There were problems both with 
the calibrator stars and with the camera's ADU conversion fac- 
tor during this night, and weather conditions were not very sta- 
ble. 



A. 2. EFOSC2 optical imaging & spectroscopy 
A.2.1. Photometry 

Broadband photometry in the B, V, and R filters was obtained 
on February 26, 2006 with EFOSC2 (Buzzoni et al. [T98"4l > at 
the 3.6m telescope in La Silla. These imaging data have a bin- 
ning of 1 (i.e. 0.157"/pix). To better resolve the emission be- 
tween the components, the data were re-sampled to 0.079"/pix. 
Eight frames of 1 sec integration time were co-added for each 
filter. Despite the short exposures, the A component was satu- 
rated in the V and R filters. The B2 source is not resolved from 
Bl in the optical filters. Components A, B, and C have a flux 
error of 0.1 mag due to the problem in resolving their fluxes. 
The aperture photometry of source C does not include the arc- 
shaped emission region. Landolt standard fields (Landolt 1992) 
for zero points were observed just prior to the science data. 
In addition, we used SExtractor to extract the photometry and 
colour determination of all ~2000 sources in the EFOSC2 field 
(about 5.5'x 5.5'), in order to check whether more objects be- 
long to the HD 155448 system. An R band EFOSC2 image is 
shown in Fig. [2] 



A. 2. 2. Spectroscopy 

EFOSC2 spectra were obtained on March 1 and April 20, 2006, 
as well as on August 19, 2007 and August 12, 2008. The spec- 
tra obtained in 2006 were taken with grism 11 and a 1.0" slit, 
covering 3380-7520 A at a dispersion of 2.04 A/pix and a mean 
resolution R~400. Slit orientations were applied as shown in 
Fig. Q] to avoid a contamination by the neighbouring compo- 
nents as much as possible and to resolve the C source and the 
arc individually. In 2007, we re-observed the C component with 
a 1.0" slit and EFOSC2 grisms 5 (5200-9350 A, 2.06 A/pix, 
R-300) and 18 (4700-6770A, 1.0 A/pix, R-600), to obtain 
higher resolution and a redder spectral coverage. The slit orien- 
tation during the 2007 campaign was east-west, to cover only 
the C component ("Slit C only" in Fig.Q}. When the high- 
resolution holographic grisms 19 (4441-51 14A, 0.34 A/pix, 
fl~3000) and 20 (6047-7147 A, 0.55 A/pix, fl~2500) became 
available, we observed the system again with a 0.5" slit in ori- 
entations as shown in Fig. Q] and Table [2] Spectrophotometric 
standard stars were taken from the list of Hamuy ( 1992 , 1994) 
and observed with a wider 5.0" slit. Bias and flats were taken 
in a standard way. The wavelength calibration was done by cor- 
relating to a HeAr lamp spectrum. The calibrated final spectra 
were averaged for consecutive exposures in each setting and 
the continuum normalised. 

Further data processing was performed for the /?~3000 spectra. 
The 3cr error in the wavelength calibration obtained is ~4 km/s 
(cr = 1.2 km/s). To produce the position velocity diagrams dis- 
cussed in Sect. 12. 31 the two-dimensional wavelength calibrated 
spectra were corrected for the curvature induced by the instru- 
ment along the spatial direction and the trace in the PSF along 
the dispersion direction. The curvature along the spatial direc- 
tion was determined by fitting a second degree polynomial to 
the centre position (determined by a Gaussian fit) of the bright 
sky [O i] emission line near 6300 A along the spatial axis. Then 
the spectra at each spatial position was shifted such that the 
centre of the sky [Oi] emission is always at the same pixel. 
The reference pixel was the pixel in which the [Oi] sky line 
crosses the peak of the star's continuum. The PSF trace was 
determined by fitting Gaussians to the PSF at several positions 
in the continuum and then fitting a second degree polynomial 
to the PSF centre positions. The spectra then were shifted such 
that the centre of the PSF is always located at the same pixel. 
The reference pixel is the median of the PSF centre positions. 



A.3. SOFI NIR imaging 

We acquired JHK broadband images with SOFI (Moorwood 
et al. [TW8j l at the La Silla NTT during May 7, 2004. SOFI's 
small field scale with 0.144"/pix was used. For each filter, ten 
exposures with random jitter box width of 20" were coad- 
ded. The total integration time in each passband amounts to 
100 sec. Darks and dome flats were taken in a standard way. 
Unfortunately, it was not possible to get separate photometry 
for the Bl and B2 components, owing to the low spatial reso- 
lution. The photometry of the C component is derived without 
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the arc, placing the aperture - as closely as possible - only 
around the star. The derived magnitudes of source C might be 
slightly fainter than they really are, since with circular aper- 
tures it cannot be avoided that the very closely located arc con- 
tributes partly to the 2 pixel thin sky annulus, thus resulting in 
a higher value for the sky. This uncertainty is considered in the 
corresponding magnitude errors. Zero points (ZPs) were de- 
rived from NICMOS photometric standard stars, which were 
observed two hours after HD 155448. We estimate a final pho- 
tometric accuracy of 0.1 mag, as we do not accurately know 
the quality of sky transparency during these observations that 
were obtained from the ESO science archive. The estimate of 
0. 1 mag is derived from the maximal difference in our ZPs to 
those of a photometric SOFI night. 

A.4. NACO NIR imaging 

Additional adaptive optics data were obtained with NAOS- 
CONICA (NACO; Lenzen et al. 120031 Rousset et al. [20031 at 
the Paranal Observatory in service mode during ESO period 75 
(i.e. between April and September 2005). We acquired images 
in narrow band filters (NB) at the central wavelengths 1 .64 /mi, 
2.12 /mi, 3.74 /mi, and 4.05 /mi, with a spatial resolution of 
0.027"/pix and total integration times up to 435 sec. Dark, flat- 
field, and sky correction were made in the standard way. Ten 
individual exposures for each filter, with a jitter box width of 
10", were shifted and co-added. 

We think that some problems must have occurred during these 
observations, because we observed strong flux losses when 
comparing the derived photometry to the expected values, as 
could already be recognised from the raw data. We immediately 
see that the exposures are not very sensitive, so that even the arc 
northeast of the C component is not seen in any of the NACO 
images. Nevertheless, we could derive useful relative photom- 
etry and astrometry for the narrow band filters NB_1.64 and 
NB_2. 12. Unfortunately, no standard star was observed for this 
service mode programme. Therefore, we first derived relative 
photometry and then scaled it to the 2MASS H and K band flux 
of the D component, since this component likely has the most 
accurate 2MASS flux, which is recognisable from the flags in 
the 2MASS catalogues, while its larger distance to the neigh- 
bouring stars will avoid contamination. The resulting NACO 
photometry matches that of SOFI within the error ranges. 

A. 5. VISIR MIR imaging & spectroscopy 
A.5.1. Photometry 

We observed the HD 155448 system between 8 and 19 /mi with 
the mid-IR instrument VISIR (Lagage et al. 120041 ) at Paranal 
in service mode in period 75. In the N band, dual imaging 
around the wavelength of PAH emission was performed with 
the VISIR filters PAH2 (1 1.25 /mi), and PAH2_ref2 (1 1.88 /mi) 
as the nearest continuum reference. The Q band exposures were 
taken with the Q2 filter (18.72 //m), which - for VISIR - has 



the best Q band sensitivity. The spatial resolution at both wave- 
lengths is 0.075"/pix in the VISIR small field. All mid-IR ex- 
posures were performed with a standard chopping and nodding 
technique with a throw of 10". Mid-IR standard stars for flux 
calibration were selected by the observatory staff from a list 
of mid-IR standard stars by Cohen (1998), and were observed 
close in time and airmass to the science observations. 

Data cube chop and nod co-addition was done in the standard 
ways. Where necessary, we corrected for known artifacts (in 
the form of stripes) of the VISIR detector. The co-added ex- 
posure time for each PAH observation amounts to 1085 sec. 
For the final image, we selected the five best out of six obser- 
vations, spread over several months, resulting in 5425 sec (1.5 
hours) total integration. The resulting PAH on/off images have 
the same exposure time and S IN. Both images were re-centred 
with an accuracy of 0. 1 pix before subtraction. In the Q band, 
each co-added observation sums to 1987 sec. For the final im- 
age, the four best observations (out of five) were selected, re- 
sulting in a total exposure time of 7948 sec (2.2 hours). 

The photometric accuracy is limited by the following uncer- 
tainties: (1) uncertainty in the aperture, which is negligible 
in our case with errors smaller than 5 mJy, (2) uncertainty in 
background subtraction: by co-addition of many exposures this 
influence got also strongly suppressed, and (3) uncertainty in 
the flux calibration. The latter may be significant, because the 
atmospheric transmission often varied in the short time be- 
tween observation of the science target and the consecutively 
observed standard star. Since several calibrated images are co- 
added, this effect may be neutralised, but nevertheless we esti- 
mate the photometric accuracy to be only within 10%. 



A. 5. 2. Spectroscopy 

Low-resolution N band spectra (R~350 at 10 /mi) were ob- 
tained between 8 and 13 /mi with a 1.0" slit. We aimed at hav- 
ing both the signal of the star and its arc-shaped circumstellar 
matter in the slit, so we rotated the slit to an orientation of 50° 
position angle (cf. Figs.|3]and[8]l. Our three VISIR spectral set- 
tings, centred at 8.8 /mi, 1 1.4 /mi, and 12.4 /urn, cover approx- 
imately the wavelength ranges 8.0-9.4 /mi and 10.4-13.1 /mi, 
after removing the outer range of each setting. 

The spectroscopic frames were debiased and checked for a cor- 
rect alignment of the spectrum with the x-axis of the frame. 
Two peaks are distinguishable in the spectral profile: the posi- 
tion of the star and that of the circumstellar arc. We extracted 
the signal separately for both regions, with "a neutral zone" left 
in between, and coadded the signal of corresponding nod posi- 
tions. In addition, data from various observations between June 
and September 2005 were averaged. An initial ad-hoc wave- 
length calibration, derived from the central wavelength of each 
spectral setting and the spectral resolution per pixel, was re- 
fined by cross-correlating with atmospheric lines. The final er- 
ror in wavelength accuracy SA/A is < 10" 3 . 
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Mid-IR standard stars for telluric correction were selected by 
the observatory staff and executed immediately before or after 
the science data. Flux was calibrated in the standard way by 
dividing the target through the calibrator spectrum and multi- 
plying the result with a profile of the standard star taken from 
Cohen (1998). Due to the varying atmospheric transparency 
in the time between observation of the science target and the 
standard star (typically less than 1 hour), the flux calibration 
is not accurate so we need to cross-correlate the spectral flux 
to our aperture photometry. Moreover, we have wavelength de- 
pendent slit losses: while we covered the entire stellar neigh- 
bourhood (silicate-dominated), the emission of the outer part 
of the circumstellar arc (PAH-dominated) was not included in 
the 1 .0" slit. Therefore, we scaled the 1 1 .4 fim silicate spec- 
trum of the stellar vicinity to the aperture photometry of 0.44 Jy 
at 11.88 yt/m (cf. Table @}, motivated by the dominance of this 
spectrum by silicates. Similarly, we scaled the 11.4 fim spec- 
trum of the PAH-dominated circumstellar arc to the photometry 
for the arc of 0.41 Jy at 11.25 yum, which was obtained as the 
difference of the total flux and the flux in the vicinity of the star 
(cf. Table 0J. We corrected most of the slit losses in the PAH 
emitting wing by this approach. Finally, the 12.4 yum settings, 
which overlap with the 11.4 /im ones, were scaled to match 
the flux of the latter. We note that this flux calibration is not a 
unique solution, but the best effort possible for this data. We 
estimate a possible deviation from the actual flux within 10%. 
More difficult is the adjustment of the 8.8 yum settings, which 
do not overlap with the more longward settings and for which 
we do not have aperture photometry at shorter wavelengths. In 
the absence of other calibration options, the 8.8 yum settings are 
adjusted to match the longward spectra. The approximate error 
in flux calibration for the 8.8 yum silicate spectrum is another 
10%, while the absolute intensity of the 8.6 yum PAH band re- 
mains uncertain. 



uum underneath the source profile. Pixels flagged by the data 
pipeline as being "bad" were replaced with a value interpolated 
from an 8-pixel perimeter surrounding the flagged pixel. In the 
case of the Short Low module, the spectra were extracted using 
a 6.0-pixel fixed-width aperture in the spatial dimension. For 
the high-resolution modules, spectral extraction was done by 
fitting the source profile with the known PSF in the spectral im- 
ages. Low-level fringing at wavelengths >20 fim was removed 
using the irsfringe package (Lahuis & Boogert 2003). The 
spectra were calibrated using a spectral response function de- 
rived from IRS spectra and MARCS stellar models for a suite 
of calibrators provided by the Spitzer Science Centre. To re- 
move any effect of pointing offsets, we matched orders based 
on the point spread function of the IRS instrument, correcting 
for possible flux losses. 



A. 6. Spitzer MIR spectroscopy 

HD 155448 was observed with the infrared spectrograph (IRS, 
Houck et al. 2004) onboard the Spitzer Space Telescope on 
March 22, 2005 (PID: 3470, PI: Jeroen Bouwman). The source 
was measured using Short Low (5.2-14.5 yum), Short High 
(9.9-19.5 yum), and Long High (18.7-37.2 yum) modules. The 
integration time was 30 seconds for the Short High module and 
six seconds for the other modules, respectively, and at least 
two cycles were used for redundancy. A high-accuracy PCRS 
peak-up was used to acquire the target in the spectrograph slit. 
The spectra are based on the droopres and rsc products, pro- 
cessed through the S18.7.0 version of the Spitzer data pipeline 
for the low- and high-resolution modules, respectively. For the 
details of the data reduction procedure we refer to Juhasz et 
al. (120101 ) and Acke et al. (120101 ). The brief description of the 
data reduction steps is as follows. For the Short Low spectrum 
the associated pairs of imaged spectra were subtracted in order 
to correct for the background emission, stray-light contamina- 
tion and anomalous dark currents. For the high-resolution spec- 
tra the background has been removed by fitting a local contin- 



